Gravitational red-shift and deflection of slow light 
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We explore the nature of the classical propagation of light through media with strong frequency- 
dependent dispersion in the presence of a gravitational field. In the weak field limit, gravity causes 
a redshift of the optical frequency, which the slow-light medium converts into a spatially-varying 
index of refraction. This results in the bending of a light ray in the medium. We further propose 
experimental techniques to amplify and detect the phenomenon using weak value measurements. 
Independent heuristic and rigorous derivations of this effect are given. 

PACS numbers; 



We live in an age where we constantly push technolog- 
ical boundaries to explore the subtleties of the quantum 
domain and to challenge our understanding of the micro- 
scopic world. Occasionally these explorations lead us full 
circle back to the macroscopic classical regime to prompt 
new and interesting questions about previously studied 
topics. These revisitations only deepen our understand- 
ing of the world. 

New materials with strong frequency-dependent dis- 
persion are also being explored that allow for remark- 
ably slow group velocities in the propagation of light 1] . 
These new developments have raised a new question: can 
we actually measure the effect of gravity on slow light us- 
ing a table-top device? The purpose of this paper is to 
describe how this could be done. 

One of the first tests of the General Theory of Rel- 
ativity was the confirmation of the bending of light 
near the sun in 1919, showing that gravity does affect 
the propagation of light 0]. However, it has been diffi- 
cult historically to test these effects in any compact way 
in a laboratory. Between the fact that gravity is a rather 
weak effect and the fact that light typically travels at 
large speeds, any measurement of a gravitational effect 
on the propagation has required a large propagation dis- 
tance - much larger than any typical laboratory. 

For Earth-confined experiments, however, the domi- 
nant measurable effect of gravity on light is the gravi- 
tational red-shift of the frequency, and not the bending 
of its trajectory (see e.g. According to this effect, 

a photon starting at an initial height and climbing to 
a final height y will have its frequency red-shifted in a 
weak gravitational field. 



uj' =uj{l ~ gy/c^). 



(1) 



The shift in frequency can be viewed as a manifestation 
of the equivalence principle: the crests of sequential light 
waves spread further apart as the beam climbs. Alterna- 
tively, the shift of the frequency can be considered to arise 
from conservation of energy as the photon climbs a po- 
tential well. As we shall see later, the height-dependence 
of the frequency continues to be the dominant effect of 
gravity during slow light propagation that also results in 
a slight bending of the trajectory. 



The shift in frequency due to gravitation may be small, 
but Pound and Rebka managed to measure and verify it 
to within 10% of the predicted value in 1959 using a 
height difference of only 22.5 meters f^. Their experi- 
ment used gamma rays emitted by the nucleus of Iron-57 
to achieve the sharp spectral lines needed to measure the 
shift. The Mossbauer effect, stating that atoms in a lat- 
tice may emit radiation from their nuclei with almost no 
recoil since the entire lattice collectively recoils, allowed 
the gamma rays to be emitted with minimal Doppler 
broadening, improving the precision. Further tests have 
since been done using greater height differences and tech- 
niques to measure the effect of gravitational redshift to 
a precision of about one part in 10^ The effect has 
even been necessary to include in the GPS guidance sys- 
tem To the best of our knowledge, the Pound- Rebka 
experiment holds the record for the shortest height dif- 
ference to see a gravitational effect on light. 

The aim of this paper is to combine the physics of slow- 
light materials with gravitational effects on the light. 
We show from fundamental principles that a light ray 
is predicted to follow a parabolic trajectory through 
the medium (to leading order in the gravitational field 
strength). The bending is predicted to arise from the 
gravitational redshift, which is amplified by the strongly 
dispersive slow-light medium. For typical experimental 
setups at the time of writing, we predict the expected 
deflection to be on the order of an Angstrom. We also 
propose an experimental set-up to detect this very small 
beam displacement. Hence, if the experiment proposed 
in this paper is realized, then the effects of the gravita- 
tional red-shift would beat the previous height record by 
many orders of magnitude. We note that this question 
has also been recently examined independently by Ku- 
mar 0. However both our qualitative explanation, as 
well as our quantitative predictions of this effect are very 
different. 



I. ESTIMATE AND INTUITION 

Later in this paper we will present a general analysis 
of the gravitational effect on light in a strongly dispersive 
medium from a General Relativity perspective. However, 
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the basic intuition of the physics can be understood quite 
simply. A slow-hght medium has an index of refraction 
that sharply varies linearly with frequency in a particular 
frequency range. A beam of light will be slightly gravi- 
tationally red- or blue-shifted as it changes height in the 
medium. The combination of these effects gives the in- 
dex of refraction an effectively linear height dependence, 
which is amplified by the steep dispersion relation. The 
height dependence of the index of refraction will then be 
translated into a spatial deflection of the light beam as 
the beam propagates. Therefore, by measuring the spa- 
tial deflection of the beam, one is indirectly measuring 
the effects of the gravitational red-shift. 

Relying on the above intuition, we can now estimate 
the size of the effect in a table-top experiment. Assuming 
that the deflection height, ?/, will be small, the index of 
refraction n{uj) may be expanded to linear order in y as 
follows. 



, , dn duj 

n(uj,y) = no + -^--j-y- 
duj dy 



(2) 



The gravitational red-shift ^ implies that dio/dy = 
—ujgjc?. The classical theory of optics [1] defines the 
phase velocity Vp and group velocity Vg of a wave packet 
as 



V 



def LO 

def duj 
' " dk 



c 
n 



. dn 

duj 



(3) 
(4) 



and using the dispersion relation k — ujn/c, we can 
c[l/vg — l/vp\. Hence, the effective in- 



rewrite uj 



diu 



dex of refraction may be written as 

/ ^ 9yc{vp-Vg) 
n{uj,y) = no ^ —■ 



(5) 



Using the approximation that Vp w c, and Vg <^ Vp, we 
can rewrite this in the simplest form that still preserves 
the characteristic velocities: 



n{uj,y) w no 



(6) 



9y_vp 

C2 Vg' 

This index of refraction leads to a deflection of the 
beam, given by solving the geometric optics equation SJ , 



d_ 

ds 



niuj,y)^ 
ds 



= Vn. 



(7) 



This is done by considering a beam initially pointing in 
the horizontal (x) direction and taking the a;-component 
of the above equation to find the conservation law, 



n((jj, y) cos — const, 



(8) 



where 6 is the downward angle of the beam from the 
horizontal. The trajectory may be found by considering 



tan 6* = — — 



dy 
dx 



' n{uj,yy 
n(w,0)2 



where the expansion is to first order in gVpz/c^Vg, and 
we have approximated no ~ 1. This simple differential 
equation can now be solved for the vertical deflection Ay, 
given a horizontal displacement L, yielding 



Ay 



(10) 



Similarly, the angle Q of the final beam relative to the 
horizontal is given by, 



gLvp 



,Ay 
' L ' 



(11) 



v„ 



The more careful analysis later in this paper confirms 
that these expressions are indeed the first order approx- 
imation to the vertical drop. 

Making an estimate, we assume plausible lab values 
of L = 0.8m, Vg = lO^m/s, and Up = 3 x 10*m/s, so 
the vertical deflection is estimated to be A « lA, and 
the angular deflection is estimated to be « 0.2nrad. 
By varying L, or Vg, the deflection can be adjusted by 
several orders of magnitude around this estimate. 



II. EXPERIMENTAL TESTS 

The nontrivial nature of this proposed experiment re- 
quires careful attention. The assumption that the in- 
dex of refraction varies as a function of height restricts 
the number of systems that can be used to deflect the 
light. For example, electromagnetically induced trans- 
parency , nonlinear magneto-optical rotations [l3| , co- 
herent population oscillations or Raman-Nath inter- 
ference in liquid crystal light valves can probably all 
achieve the experimental products needed for an observ- 
able signal. However, all of these systems use a strong 
pump beam to prepare a coherence or interference in the 
medium. The pump beam is also experiencing the grav- 
itational redshift. For pump beams of nearly the same 
frequency as the probe, the redshift of the probe is then 
washed out, because the index of refraction established 
by the pump beam also shifts in the same direction as 
the probe beam. Since dyUj oc uj, only when the pump 
beam is dramatically different in frequency will a pump- 
prepared system have the ability to detect the red-shift. 
In addition to the problems associated with a shifting dis- 
persion, pump beams with a nonuniform intensity profile 
(e.g., Gaussian beams) also cause some guiding effects 
due to index of refraction changes. Some absolute fre- 
quency systems also have technical difficulties. For ex- 
ample, hot atomic vapors have gravitational and thermal 
density gradients that would usually overwhelm the small 
deflection [13,]. 

One can envision certain classes of systems that have 
promise to measure the deflection. One class of systems is 
one in which a pump beam and probe beam propagate in 
opposite directions. This is effectively the same as having 
a pump beam blue shifted while the probe is red shifted. 
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The deflection would be twice as large, because the dis- 
persion is shifting in the opposite direction in frequency 
space as the probe as a function of height. Another class 
of systems is one which uses fixed frequency, narrow res- 
onance, solid state, bulk materials. If the inhomogeneous 
broadening is negligible, the probe frequency can be very 
close to the resonance frequency to achieve small group 
velocities. These experimental parameters are demand- 
ing. 

In order to measure these small gravitational effects, 
we enlist quantum mechanics for help. Interference ef- 
fects can measure length differences of much less than 
a wavelength. Indeed, Manly and Page have proposed 
to essentially directly measure the change in the index 
of refraction ^ by interfering a red-shift beam with a 
blue-shifted beam in a fiber-based Sagnac interferometer 
[l3 |. Such an approach can also be taken in a slow-light 
medium. We describe below another quantum approach 
using another kind of quantum interference. In order 
to measure such a small defiection of a light beam, we 
propose to extend recent advances in precision metrol- 
ogy by amplifying the deflection using quantum "weak 
value" measurements [l3|- Weak values have the prop- 
erty that the post-selected expectation value of an op- 
erator can exceed the eigenvalue bounds of the operator 
[igL 113] ■ Recently, Hosten and Kwiat successfully used 
optical weak value measurements to detect a 1 A deflec- 
tion in an optical beam ^18]. In the Hosten-Kwiat exper- 
iment, the weak value operator was the polarization of 
the light beam, which was entangled with the transverse 
position degree of freedom. Most slow-light materials do 
not have a polarization-dependent index of refraction, so 
we propose to use instead a 'which-path' operator of an 
optical Sagnac interferometer 0] as shown in Fig. 1. A 
laser source is at the input of a Sagnac interferometer. If 
no other optical element is in the interferometer and with 
ideal alignment, all of the light will exit the input port 
of the interferometer. The optical path length of both 
directions are identical in a Sagnac because both paths 
simply traverse the same route but in different directions. 
The reason all of the light exits only the input (or bright) 
port is due to a relative 7r/2 phase shift for each reflection 
versus transmission through the beam splitter, resulting 
in exactly destructive interference for the dark port. Us- 
ing a tunable birefringent element (a Soleil Babinct com- 
pensator, SBC), we can break the symmetry and cause 
a relative tunable phase between the two directions in 
the interferometer. For this gedanken experiment, we as- 
sume the light entering the interferometer is horizontally 
polarized. The light that is propagating in the counter 
clockwise direction is rotated to be vertical via a half 
wave plate oriented at 45° with respect to the horizontal 
polarization. We align the SBC such that the vertically 
polarized light receives a relative phase shift compared 
with the horizontally polarized light which is propagating 
in the clockwise direction. The light propagating in the 
clockwise direction is then made to be vertically polar- 
ized. The two paths then interfere at the beam splitter 



Dark port 



50/50 



LASER 



Waveplate 




Slow Light 
Medium 



SBC 

FIG. 1: Experimental Setup. A laser beam enters a Sagnac 
interferometer. If no other optical element is in the system 
and for perfect alignment all of the light exits out the input 
port of the interferometer and goes back to the laser. Using 
a combination of a half (A/2) wave-plate and a Soleil Babinet 
compensator (SBC), one can tune the relative phase acquired 
for each direction in the interferometer. This allows us to 
tune between the dark and bright ports, which is necessary 
for weak value measurements. The slow light medium the 
interferometer allows us to couple the 'which-way' informa- 
tion of the interferometer to the transverse deflection of the 
light beam. A quadrant detector, which measures beam de- 
flections, is in the dark port of the interferometer. 



but the counter clockwise light has picked up a tunable 
relative phase shift of (j>. 

A slow light medium (SLM) is placed in the interfer- 
ometer, causing a transverse gravitational deflection Sy 
of the beam in the downward direction. We want to 
amplify this deflection using weak value techniques. By 
placing the SLM at an asymmetric point in the optical 
path, we can arrange that there is a long path of length 
Li (corresponding to the clockwise direction O having a 
larger deflection Syi) and a short path of length Lg (cor- 
responding to the counter-clockwise direction O having a 
smaller deflection Sys) between the SLM and the 50/50 
beam splitter. We now describe how to use weak val- 
ues to amplify the small spatial deviation. We define 
the system as being the which-path information of the 
light beam, and the meter as the transverse profile of the 
beam. The system operator is the which-path operator 
A = (5y;| 0)(0 \ + Sys \ 0)(0 |. Following the propagation 
of the state through the Sagnac interferometer, starting 
from a pre-selected system state \tpi) (formed at the BS), 
and postselecting on a system state | ■;/'/) (exiting the dark 
port) we find that the weak value 



(^/|A|V'.) 



of A is given by a purely imaginary result, 
. _ . 5yi - Sys 

J±qn ~ I : . 



(12) 



(13) 
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The SBC phase <C 1 characterizes the deviation from 
perfect darkness, and is res pon sible for the amphfication 
of the deflection. See Ref. [l^ for detailed calculations. 
Expressing the deflection in terms of the gravitational 
deflection angle 9 (fTTj) . 5yi — 5ys = [Li — Ls)9, we have 
for the total deflection, 



(Sy) 



F{Li - 



(14) 



where is a factor originating from the dynamical prop- 
agation, which brin gs a dditional enhancement (a factor 
of about 100 in Ref. [18|). This deflection is then directly 
measured with a quadrant detector in the dark port. 



III. HAMILTON'S OPTICS 

In order to give a rigorous treatment of light prop- 
agation through media with anomalous dispersion in a 
gravitational field, we start with the wave equation it- 
self. Many of the historical treatments of light propaga- 
tion have implicitly assumed the index of refraction to 
be weakly dependent (or not dependent) on frequency 
since the derivations were done before slow light was a 
commonly studied phenomenon, e.g. Ref. [23|. Hence, a 
revisitation of the physical foundations is illuminating. 

The components of an vector EM wave propagating 
through a material medium obeys the modified wave 
equation. 



(15) 



Here 5* is any component of a vector EM wave, n is the 
index of refraction of the material, and c is the speed of 
light in vacuum. 

We assume the propagating wave components have the 
form ^ = Aexp{iS/X), where S{f,t) is the eikonal of 
the wave [81 . Then, to leading order in 1/A, we get an 
equation for the eikonal itself in the geometrical ray ap- 
proximation. 



dt 



\ys\' 



0. 



(16) 



The eikonal relation holds for all components of the vec- 
tor EM wave so the vector nature of the EM wave itself 
is unchanged during the propagation as a ray. By taking 
the leading order in 1/A, we assume A is small compared 
to the length scales through which the wave propagates, 
such that the wave behaves effectively as a ray. Note 
that the ray may still undergo small angular deflections 
without violating this approximation, which is a crucial 
point for the results in this paper. 

In small regions of space and time, the eikonal can be 
expanded in a series to first order in space and time. 



S So + f-VS + t 



as 



(17) 



where we have incorporated the wavelength A into the 
eikonal for notational simplicity. Comparing this approx- 
imation to the eikonal for a plane wave, S ^ k- r — Lut + a, 
we define the frequency and wave vector of the wave for 
a local region where the eikonal can be approximated as 
planar, 



k VS, 



OJ 



def 



dS_ 
df 



(18) 
(19) 



We then rewrite the eikonal equation ([16)) in terms of 
the local wave vector and frequency for each infinitesimal 
region of space and time. 



k-k - (^^yn^{uj,f) = 0. 



(20) 



Here the potential dependence of the index of refraction, 
n, on the local frequency and position is written explicitly 
as a reminder. 

The eikonal equation can be rearranged and square- 
rooted, giving the usual dispersion relation for the 
medium in a localized region. 



c\k\ = ujn(LU, r). 



(21) 



The signs are chosen to keep the frequency, uj, positive. 

An analogy can be made between the mechanics of 
material particles and the motion of wave packets that 
exposes a new set of canonical variables and allows a 
systematic derivation of the equations of motion for the 
wave packet. 

From Eqns. [T^ . we see that the eikonal, S, plays 
the role of the action of the system, S, with the wave 
vector, k, corresponding to the momentum of the wave, 
p, and the frequency, uj, corresponding to the energy, H. 
In particular, as is well known from Quantum Mechanics, 
the energy and momenta of a wave packet composed of 
plane waves are directly proportional to the frequency 
and wave numbers, respectively. The scaling constant 
in Quantum Mechanics is H, implying that the classical 
action is just the scaled eikonal for the wave packet. 

The explicit correspondence for the translation to 
Hamiltonian dynamics is 



S ^ S 

p ^ k 

t ^ t 

H ^ LO. 



(22) 



Using this correspondence, either the eikonal equation 
(PO)) or the dispersion relation (|2ip can be understood to 
represent Hamilton-Jacobi equations of the system in an 
implicit form |2l| . Solving each for lo would give the stan- 
dard form of the Hamilton-Jacobi equation using these 
identified canonical coordinates. Furthermore, w will be 
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a conserved quantity since the equations are independent 
of t. For convenience in referring to these equations, we 
wiU denote the eikonal relation as F and the dispersion 
relation as G: 



F{k,uj) =^ k-k 



G{k,uj) =^ \k\--n{uj,r)^0. 
c 



{^^)\\Lo,r)^Q, (23) 
(24) 



Using the conserved quantity uj as the Hamiltonian of 
the system, we can write down Hamilton's equations of 
motion directly in the form. 



dk ^ 
—r — k 
dt 



V, 



clef dr 
' " di 



dk' 



(25) 
(26) 



It becomes clear at this point that Vg must be interpreted 
as the group velocity of a wave packet that obeys the 
eikonal equation, and not the phase velocity of the car- 
rier wave itself. Hamilton's equations of motion along 
with the eikonal equation completely determine the mo- 
tion of a wave packet treated as a point-particle through 
a material medium along a geometrical ray. The geo- 
metrical ray approximation made to derive the eikonal 
relation implicitly assumes wave packet behavior. 

If the form of the index of refraction n is known, then 
F or G may be solved for uj, and Hamilton's equations 
can be solved for the trajectory of the packet directly. 



IV. INFLUENCE OF GRAVITY 

In order to include the influence of gravity on the slow 
light beam, we recall that the main gravitational effect 
for earth-based experiments is the gravitational red-shift. 
The following simple rule can be used to make a weak field 
translation of the dispersion relation to include gravita- 
tional effects. 



(27) 



where U = —GM/c^r is the scaled Newtonian gravita- 
tional potential. We give a general relativistic justifica- 
tion of this rule in the appendix. It is the spatial depen- 
dence of U in the shift of the frequency that will cause a 
bending of the trajectory of light. 

The gravitational redshift is very small, ?7 <C 1, so we 
expand the dispersion relation after the replacement 
([27)1 to first order in U, 



G(k,r) 



« lfcl 


Lon 




n + UJ 




c 


c 




« lfcl 


--( 






Vp \ 




= 0, 









dn{Ld) 
duj 



(28) 



where we have used the standard definitions of Vg and Vp 
pi4p . We will see in the next section that the definition 
of Vg is effectively unchanged. This final simplification 
hides the simple frequency redshift that is occurring, but 
makes the form of the first order correction in terms of 
the quantity Uvp/vg particularly apparent. 



A. Hamilton's Equations 

We proceed to solve Hamilton's equations of motion 
with the frequency uj playing the role of the Hamilto- 
nian of the system, and the new dispersion equation (|28p 
acting as the Hamilton- Jacobi equation. 

Hamilton's equation (|26|) gives the vector group ve- 
locity of the packet. The frequency uj involves only the 
quantity |fc|, so ^ oc fc. Therefore, we now focus on the 
magnitude of the group velocity, 
derivative to (^5)1 yields 



Applying a frequency 



d\k\ 
duj 



{l-U) 



n{uj) 



duj 



U 



dn{uj) 

duj 
ijj^ d'^n{uj) 
c duj'^ 



(29) 



Now we make the slow light approximation that n is 
linear in uj. Slow light occurs for a frequency band where 
the index of refraction n has a steep linear dependence 
on UJ with positive slope. Therefore we neglect the last 
(second-order derivative) term in Eq. (|29p . so 



df 
'dt 




(1 - U)n + uj{l 



(30) 



We see that the group velocity is modified very slightly 
from its traditional form by the presence of the gravita- 
tional field. However, the packet will always travel in the 
same direction as the wave vector. 

Hamilton's other equation, (j25p . will determine how 
the wave vector itself changes over time and thus how 
the trajectory will curve. Since the only position depen- 
dence in w is radial due to the presence of U, the gradient 
reduces to a radial derivative. We can implicitly find drUj 
by applying a spatial derivative to G, 



d_ 

dr 



c\k\ = {I - U)Lon - uj'^U^ 

auj 

o.|i(a-.,„.a-3.,.£ 



dr 



- UJ 



dn 
duj 



Solving for drUj gives 

duj 
dr 



dn 



(1 - U)n+{1 



3C/)^iS 



dU 
dr 



(31) 



(32) 
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Then from p5|) we find the change in wave vector to be 



dk 



-Vlo 



duj 
dr 



. dn 



(f - U)n+il 



duj 



Now we linearize the dimensionless gravitational po- 
tential U to give a more useful form for a lab, and eval- 
uate the derivative in (1331). 



U = 



MG 



-gR 



'■{R + r) 



so, 



dn 
duj 



(l-C/)n+(l-3t/V^ 
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9 - 



(34) 



(35) 



Thus we see a net deflection of the wave vector having 
to do with the presence of the local gravitational accel- 
eration g. Here R is the radius of the Earth, M is the 
mass of the earth, G is Newton's gravitational constant, 
and r is the height displacement from the initial starting 
position of the packet. Notice that the scaling factor in 
([55)1 is the ratio of the new group velocity to the normal 
group velocity. 



Noting the fact that g/c 



10 



gRic 



10^^^, and ijj ~ 10^°s^^ for the surface of the Earth and 
the optical domains we are considering, we can reduce 
equations ((30)) and ([35]) to simpler final forms. 



k w ttIj 



(36) 



(37) 



The group velocity appears unchanged to a very good 
approximation with such a weak field from what is tra- 
ditionally expected in a slow light medium, but the wave 
vector picks up a small, effectively constant, deflection 
over time that is dependent on the carrier frequency of 
the packet and the local gravitational acceleration. The 
deflection points toward the gravitating body, so light 
will bend toward the body, matching intuition. 



B. Wavepacket trajectory 

With the newly derived approximate equations of mo- 
tion for the weak field limit at the surface of the Earth, 
([5S)) and (133), we can solve the trajectory of a point-like 
wave packet in a straight forward manner. 

Assume that for a small local bit of trajectory the ra- 
dial direction can be well approximated by the vertical 
Cartesian direction y, and consider the wave vector point- 
ing entirely in the horizontal Cartesian direction x at an 
initial time ii = 0, so. 



MO) 



-X. 



(38) 



Here Vp is the phase velocity of the packet. 
Solving for k{t) from ([571) yields 



k{t) 



^9.. 

—ty 



m 

9 



ty 



(39) 



so 



kit) 



k 



1/2 



9 

X - Vp—ty. 



(40) 



Thus the wave vector rotates over time, and its magni- 
tude increases due to the gravitational redshift. Notice 
that the wave rotation is not dependent on the group ve- 
locity, but only on the carrier phase velocity. The group 
velocity only appears when we calculate the position of 
the packet. The final approximation (|40|) is made since 
we are only concerned with the lowest-order drop [2^ . 

Solving for r(t) from ([M)) . starting the packet at an 
initial position r (0) at the start of a medium of length L 
yields 



k{t)dt, 



f(0) + Vgtx - ^^VpVgt'^y, 



(41) 



which is a parabolic trajectory analogous to a kinematic 
trajectory, but with an effective gravitational accelera- 
tion given by. 



def VpVg 
9eS = 5-V 



(42) 



Thus a wave packet drops much less due to gravity in a 
medium with slow group velocity than it would in vac- 
uum in the same amount of time, explaining why we have 
not casually observed them dropping like rocks in previ- 
ous experiments! 

Since we are interested mostly in how the wave packet 
trajectory is altered after traversal though a finite slow- 
light region, we set the final horizontal position to the 
length of the propagation region: rx{tf) — L. Solving 
for tf, we see that 



L 

tf - — ■ 

Vn 



(43) 



We can now solve for the vertical drop after a traversal 
of horizontal distance L, yielding. 



A 9 VpL^ 



(44) 
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Again, this drop is analogous to the drop found in a kine- 
matic trajectory, but with the effective gravitational ac- 
celeration (|42|) and horizontal velocity Vg. Despite the 
fact that the effective gravitational acceleration is re- 
duced, the long propagation time compensates for this, 
resulting in an amplified deflection. 



V. CONCLUSION 

We have shown that a slow light medium is expected to 
amplify the effects of gravity on a wave packet such that 
they may become visible to a sensitive apparatus using 
weak value measurement. The primary weak field gravi- 
tational effect is the induced position-dependence of the 
carrier frequency of the light, uj —^ (1 — U)uj, where U 
is the dimensionless Newtonian gravitational potential. 
The position-dependence of the frequency shift is ampli- 
fied by the strong frequency dependence of the slow light 
medium, which translates to a magnified bending of the 
trajectory of the light as it moves through the medium. 

Using the modified dispersion relation for the EM wave 
components of the light ([28|l as a Hamilton- Jacobi equa- 
tion, we applied Hamilton's equations of motion to find 
a parabolic trajectory of the slow-light pulse through the 
medium ((4T|) . We found that the trajectory is analo- 
gous to a simple kinematic trajectory, but with a much 
smaller effective gravitational acceleration. The smaller 
effective acceleration is more than compensated by the 
longer travel time, leading to the amplified vertical drop 
(j44|). For typical current lab values, we expect a vertical 
drop on the order of f Angstrom for a traversal distance 
of roughly Im. 

Although the predicted displacement is tiny, it bor- 
ders what is currently possible to detect with the am- 
plified precision of a weak value enhanced measurement. 
We have proposed an optical weak values measurement, 
based on post-selected weak measurements of a which- 
path operator in a Sagnac interferometer. This proposal 
has the advantage that in contract to previous experi- 
ments [ill, the optical deflection does not need to be 
polarization dependent. We stress that this fact makes 
the idea applicable to amplify an optical deflection orig- 
inating from any optical element. Aside from the fun- 
damental aspects of detecting a gravitational red-shift in 
a laboratory setting, a table-top detector would open a 
new field in precision gravitational metrology. 



APPENDIX A 

The purpose of this appendix is to describe how to 
relativize slow-light equations of motions. We will con- 
vert the Eikonal equation into relativistic form, then pass 
to curved space-time, and then specialize to the weak- 
field case in order to justify the red-shift replacement 
a; — > (1 — U)uj. We will also comment on other methods 



that we have used to independently check the results in 
the paper. 

Flat Spacetime. — The equations developed for Eu- 
clidean flat space parametrized by a universal time can 
be generalized to a Minkowski spacetime in a straight- 
forward manner. The generalization involves moving 
from the flat product space K x with Euclidean metric 
on each piece to the flat spacetime manifold M with a 
pseudometric having signature (— , +, +, +). 

To convert the flat space eikonal relation (PO)) into the 
language of spacetime, we note that the eikonal is a scalar 
field on the manifold. The natural replacement for the 
frequency and wave number of the wave is the oneform 
field (written in boldface) defined as the exterior deriva- 
tive of the eikonal, 

k dS = — dx^ = fcudx^. (Al) 

We use Einstein summation convention so all doubled 
Greek indices have implied sums over all 4 dimensions 
of spacetime, and Latin indices have implied sums over 
only the 3 spatial dimensions of spacetime. 

Written in the component representation with the co- 
ordinate basis the wave number oneform field is 

f dS dS dS dS\ _ f _uj , , , \ 

Here the notation = indicates an isomorphism to 
through a particular component representation for con- 
venience. In order to properly translate the eikonal equa- 
tion (j20p to spacetime, we will treat the appropriate com- 
ponents as the components of this oneform field. 

Keeping in mind that the dot product in the eikonal 
equation (|20p should correspond to application of the 
metric, we get the following translation into spacetime, 
in component form, 

F{k, r) rj'^hkj + ?7°°fcofcon2(fco, r) = 0. (A3) 

Notice that this equation is not Lorentz invariant un- 
less n = I. The presence of the medium breaks the sym- 
metry and forces a preferred frame for interpretation of 
the components of the wave vector. This is precisely 
why the components have to be interpreted in the cor- 
rect frame for the lab measurement. Assuming the lab 
is stationary, the correct frame is the coordinate frame. 
Changing to different frames results in frequency contrac- 
tion (Doppler) effects from the lack of Lorentz invariance. 

Curved Spacetime. — The transition to curved space- 
time with metric g follows by considering the spacetime 
version of the eikonal equation (|A3p to remain valid for 
any orthonormal frame and co- frame ^23.] . It is always 
possible in a curved spacetime to find an orthonormal 
frame and co-frame, though they will not be a coordi- 
nate frame pair. 

The key assertion is that the spacetime eikonal equa- 
tion (jA3p is valid in this new orthonormal co-frame 
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{kn k^). Finding the proper form of the equation 
for the lab frame requires performing a frame change. It 
is important to interpret the components of the one form 
field k as a measurable frequency and wave number only 
in the correct frame-namely the lab frame-due to the fact 
that the components will change with the frame. Note in 
particular that there is no coordinate change happening 
during this process. Everything on the manifold will still 
be parametrized by the original Cartesian coordinates. 

The frame-change results in a frame-transformation 
matrix, e, that is related to the Minkowski metric, 77 , 
the metric of curved space, and the oneform field, k, 
as, 

(A4) 

(A5) 
(A6) 
(A7) 

If the form of the metric is known for a particular lab co- 
frame (usually the coordinate co- frame), then the form 
of the transformation matrix e is completely specified 

Using these the above relations to transform the 
eikonal equation to the coordinate co-frame, we find that, 

F(k,f) r^'^kk+v"%~kon^iko,r^, 



kfi 






= {e-'Yk-^, 




— p'^p'^n 







{{e-'te^o) {(e-XO ^'(egfcp,r)] 



0. 



(A8) 



The wave oneform components in the lab frame are now 
precisely in the form we can interpret physically as in 
Eq. (jA2p . Note that the matrices from the co-frame 
changes do not completely cancel due to the lack of 
Lorentz invariance in the equation. Furthermore, the in- 
dex of refraction gains a spatially-dependent part. All 
the components of the wave oneform become mixed in 
the equation on frame changes, leading to various fre- 
quency contraction (Doppler) effects in different frames. 

In the special case when n = 1, the co-frame change 
matrices do exactly cancel and the eikonal equation be- 
comes Lorentz invariant. This special case is the tradi- 
tional case for vacuum propagation, Fn=i — k^k^ = 0. 

The new eikonal equation (jA8[) gives us a Hamilton- 
Jacobi equation for the propagation of a wave packet 
through curved spacetime in the geometrical optics ap- 
proximation in the lab frame, assuming the medium 
through which it propagates has no other effect than to 
introduce the index of refraction into the wave equation. 
In particular, any microscopic interaction in the medium 
which would prevent the packet from falling due to grav- 
ity is being ignored. 

As before, if n is independent of t, then the entire 
eikonal equation is independent of t. Therefore, w be- 
comes a conserved scalar quantity of the motion. It is 
known as the "world frequency" since it is invariant with 



respect to the world time. The "proper frequency" is 
conjugate to the proper time of an observer, rather than 
the world time, and differs by a frame transformation 
matrix from the coordinate frame to the observer frame. 

Weak Field Limit. — In order to compute some- 
thing measurable, we use the weak-field limit of the 
Schwartzchild metric for a radially symmetric gravitating 
body (the Earth). In Cartesian coordinates the non-flat 
part of the weak-field metric and the transformation ma- 
trix from flat spacetime are given by the relations. 



ief 2GM 



9m 

„-i 



2r' 
-{1-Uf 
-{l + Uf 
diag(l + [/,!, 1,1), 
diag(l-C/, 1,1,1), 



-(1-2(7), 
-(1 + 2(7), 



(A9) 

(AlO) 

(All) 
(A12) 
(A13) 
(A14) 



where Vg is the Schwartzchild radius of the gravitating 
body, and U is the dimensionless Newtonian gravitational 
potential. The rest of the metric g is identical to the flat 
space metric "q in this weak-field limit. 

The eikonal equation (|A8[) then reduces to the form: 

F(k,f) .9"^fc,fc4((e-i)j,er)((e-i)^^e^^) + 
{{e-')le'',){{e-^)Y,) n^e^^k.^r)] 
= .g'^fc^fc, +.9""fcofco"'(e[!fco,^ 
= fc.fc-(l-2C/)(-^^)'n2(-(l-C/)^,f) 



^^""'^ n2((l-[/)c.,r) (A15) 



k-k 




The last step simply omits constants in the functional 
form of n, showing the relevant dependence, and puts 
the gravitational dependence in a more consistent form. 
Notice that the simple diagonal form neatly keeps the 
frequency components unmixed, and only leaves gravita- 
tional influence attached to w as a small correction. Since 
the dimensionless potential U is dependent purely on the 
radius, this equation is invariant in form with respect to 
a coordinate change to polar coordinates as well. Tak- 
ing the square-root of this equation gives the modifled 
dispersion relation, justifying the gravitational red-shift 
replacement ([Tf)). 

We now briefly describe other independent checks done 
on our results. Rather than solve Hamilton's equations 
in three dimensions, it is also possible to give a four di- 
mensional derivation. In the 4D approach the frequency 
is conserved, and that quantity F or the quantity G are 
treated as the Hamiltonian and an arbitrary global pa- 
rameter r is introduced. The equations of motion even- 
tually reduce to the 3D ones ^2^. 
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